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I , Introduction 

This Memo describes some of the more important properties 
and manipulations of Hypergeometric Functions, P (.1, which 

p q 

may be useful in MATHLAB. A convention for representing 
F {.) is adopted which is readily adaptable to LISP operations. 

The most general type of HGF with which we will be con- 
cerned is a function of a single independent variable, x, and 
is parametricized by an "A M list, of length p, and a "B" list, 
of length q. The latter consist, in general, of atoms; the 
argument is usually x, but may also be a simple function of x. 

The power of the Hypergeometric Formulation arises due to: 

(i) Nearly every Tabulated Function can be represented 
as a HGF, for example the entire families of Bessel functions, 
and Orthogonal Polynomials, as well as Elliptic, Error, and 
Incomplete Gamma Functions* 

(ii) Linear operations, such as differentiation, integration 
(definite or indefinte), integral transforms, etc*, amount to 
simple manipulations of the list representations of the HGF. 

(iii) Analytic continuation is easily accomplished by opera- 
tions which allow the variable, x, to map into 1/x, x/(x-l), 
1-1/x, 1/(1— x) , simultaneously with simple linear operations on 
the "A" and rt B" lists. (This statement should be qualified to 



the extent that such continuations are a lwa ys possible for 
2 F 1 (-), are not required for p F-(0 where both p and q are 
less than two, and are not known — possibly do not exist — 

when both p and q are greater than unity. 

(iv) The solution of a wide class of homogeneous DE's 

with nonconstant coef ficients, say, of the form: 

(ax 2 + bx + c)y" + (dx + e)y' + f y = 

can be represented as a HGF, after an appropriate change of 
variables. 

Other advantages of the approach could be enumerated, as 
should be clear. In essence, this is also the principal draw- 
back of the technique so far as user-oriented applications are 
concerned* 

Itemizing: 

(i) It is time* and core-consuming to recognize and interpret 
properly into HGF-form an arbitrary function presented by the 
user. Once this has been accomplished the lowest-level operation, 
are rapid, as is claimed above. However, the result of the latte: 
manipulations muse also be represented in HGF notation, and now 
the process has to be "unwound". In general, this amounts to a 
fairly long sequence of tests for various linear relationships 
among members of the "A" and "B" lists, which, if satisfied, 
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identify the HGF as a known tabulated function. If these tests 
fail, then the HGF must be returned as is, {It snould be men- 
tioned that if this happens to be the case, the user now has an 
easily computable series representation,) 

(ii) In the same vein, such an all-or-nothing approach may 
be foolhardy. For example, if we wished to evaluate: 



;• 



xe x dx ( 

we could represent"~it as 



s 



- jF f-l;arax)) *( p F (a,b,c, . ..;a,b,c, . ..?x))dx . 

It works, but who needs it? 

(iii) Algebraic and rational functions are difficult and 
awkard to represent in terms of D F a (-)- Thus, when these are 
present in a function presented by a user, they must be sorted 
out and dealt with in a different manner* 

(iv) Similarly, elementary transcendental functions (i.e. 
log(.)t and all of those based on e ) are not to be subject to 
transformation into HGF representation. 

(v) It seems from the above that the user must be subject 
to some fairly strict, though not overly demanding, ground rules 
as regards the convention of the ordering of the elements of the 



, 3 - 



function presented- Ail we ask is that these elements be 
ordered according to the following hierarchy; 

(algebraic, rational) (log, e x , sin, cos, etc.,} (higher fns.) 
For example, if the given function were 

F(x) ■ f^x) + f 2 (x) + ... + f n tx) < 

then each of the fv^ x ) would be represented by an ordered triple, 
e.g. ■ . 

F{x) » x k sin x I {x) + log(x) 
((x**k)(sin x)(£hyp (1) ( (x**2)/4) ) ) 
+ (O(log x)()), 

where I Q (x) is the modified Bessel Function: 

V x) " o P l ( rl?x2/4 >' 

This format, or something similar (I am open to suggestions}, 
should assist the interpreter considerably and help to avoid 
errors. 

(v} concerning Differential Equations {see I (iv} above) , 
the fundamental solutions are indeed linear combinations of 
-il'it-J's: there are other solutions, of course, depending on the 
solutions of the associated indicial equation. However, these 
can be constructed in a straight-forward manner by means of the 
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transformation and continuation formulas mentioned in X (ii) 
and I (iii) above- similar comments hold for the Hypergeometric 
DE which yields j^f.) as its fundamental solution. For other 
values of p and q ( the question is usually rarely encountered, 
trivial, or unanswerable. 

It should be made clear that the purpose of the research 
effort outlined herein is not intended to supplant or to rewrite 
substantially the present MATHLAB; rather, the aim is to augment 
the present capability for symbolic manipulation in order to deal 
with a wider class of functions which are often encountered in 
applications, and which are not recognizable by the present 
system. Indeed, as implied in the comments under "drawbacks", 
above, such methods should be employed only as a last resort, 
when other more conventional techniques have fallen short. It is 
intended to be another tool in the kit, as it were, rather than the 
be-all and end-all. In its defense, however, it may be noted 
that the example given in "drawbacks" (v) above* can be dealt 
with efficiently only be recourse to the HGF notation and opera- 
tions, if any claim to generality is to be realized. 



♦This expression, or its equivalent, occurs with appalling fri 
quency in applications as diverse as MHD, Quantum Mechanics, 
and statistical Detection Theory. 
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The Theory of Hypergeometric Functions is quite extensive, 
and on the surface, at least, complex. Yet at the same time, 
it is very compact and very nearly trivial once all of the mathe- 
matical superstructure, accoutrements, and other random "ele- 
gances' 1 and similar rubbish have been cleared aside. Essentially, 
an understanding of the entire body of proofs requires no more 
than a rudimentary familiarity with Taylor's Series and the con* 
vergence thereof, and the technique of Analytic Continuation. 
The rest is Algebra. The most difficult and subtle task is to 
establish a basic set of "axioms 1 ' from which all other properties 

and operations may be derived. The search is still going on, but 
it appears that they are at least four in number, and perhaps as 
many as seven. 

The principal references used below are Slater (1), slater (2), 
Bateman {3), and Bureau of Standards (4). 

_II» Definition 

The Hypergeometric Function is defined by the Series: 



pWV*3 V b l ,b 2' b 3 V X) 

- F (A;B;X) 

p q 
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« (. 1 , a (. a ) n (.... )<. p ) n „* 
" ) — « (1! 

where 

(di ,a 3 , .. -,a ) ■ A, is the "A 1 * list, of length p, 
(h 1# b 2 ,...,b ) « B, is the "B" list, of length q, 
and 

(c) k - (c+k-1) V(c-l): 

= (c)(c+l)(c+2).„(c+k^l) > 
where we assume that k is an integer, and x is the indepen- 
dent variable* where c is not an integer, the second defini- 
tion of (c) v can always be used; the first definition can be 
used if the factorial functions are replaced by the Gamma func- 
tions of (c+k) and (c) , respectively . In various applications, 
it may be of advantage to adopt one of the three. 



III. Elementary Properties 

(It should be pointed out beforehand that in nearly all 
"interesting cases*' either p=q or p=q+i.) 

( i) Convergence : 

(a) if p = q, the series converges uniformly for all 
real x. 
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(b) if p = q + 1, the radius of convergence of the 
series is jx|< If however, analytic continuation 
can generally be invoked to include the region 

xj > 1? the singular point lies at x = 1; thus 
either series converges for lx| =1, arg(x) / 0. 

(c) if p > q + l t the series does not converge in any 
usual sense*. 

(d) if p < q, the series converges uniformly for all x. 
(ii) Trivial Reductions: 

(a) if one or more members of the "A" list are zero, 
the value of the series is unity. 

(b) if one or more members of the "A" list is a negative 
integer, say -N, then the series reduces to a finite 
polynomial of degree N. 

(c) if any member of the "A" list is identical to any 
member of the "B" list then both of these are 
deleted from their respective lists, and both p and 
q are reduced by unity, e.g. 



♦There are ways of handling this situation, however, by recourse 
to the theory of MacRobert's E-function or Meijer's G-function; 
we will not consider these at present. 
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2 F 2 (a,b;b,c r x) «= ^{ajcrx) , 
3F2(1.2,3r2,3;x} - ^(1, ;x) . 

(d) eyerv. iFo (.) is a simple algebraic function: 

CO 
!F <a; ;x) = £ ^ X* 



n=o 
d-x)" a . 



(2) 



(e) every ^ ( . ) i s a Bessel Function: 

oV '*'*> "WW***). 



x 2 



OV -t-l =% iJa .!(2xS. 



(3) 



x 2 



(f) if one or more members of the "B" list is zero or 
a negative integer, and cannot be cancelled by an 
an identical member from the "A" list, as in (c) 
above, then the series is undefined. 

It is worthwhile to digress briefly in order to amplify some 
of the consequences of this restriction {f ) . it must be empha- 
sized that when at least one of the b^ is zero or a negative inte- 
ger, we cannot represent the function in power-series form mean- 
ingfully - no more, no less. This may or may no t have a bearing 
on the question as to whether or not the function exists in some 
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other sense. The purpose of this aside is to illustrate some of 
the subtleties of the decision processes involved in the HGF 
formulation, and to point out that it has some limitations of an 
unexpected kind. 

Consider the function P(s;n) defined by the integral: 

(which happens to be the Characteristic Function of Fisher's 
F-Distribution (4) J • There is no doubt that the integral converges 
for n^i and Re(s)>0, but does it possess a Taylor's Series? 
Well, sort of. 

Let us attempt to formally expand F(s?n) as 



(j^) F(s;n) 



k: 



t-0 

It is not difficult to shew that this turns out to be 

FlBl * V (n-l)(n-2)...(n-k) k: 

m 7"" (n) (n+1) . ..(n+k-U s^ 

hg* (-n+1) (-n+l+l)...(-n+k) k: 

" ^ (l-n) k kl 

= -F.lnrl-nrs), by "definition", at least in 

the purely formal sense. However, as is easily seen from the 
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first line above, the term in the series for k ■ n blows up, as 
do all subsequent terms. The whole point of this exercise is to 
show that a HGP with one or more zero or negative integer members 
in its "B" list may or may not represent a meaningful function. 
To repeat, this simply means that no valid Taylor's {or Laurent) 
series exists. However, integral representations often exist in 
such cases, as the present example demonstrates. A subsequent sec- 
tion will list some of the more important of these. 

At the risk of appearing to lapse into pedantry, one further 
comment must be made — and this is probably the most important of 
all. Most of us with a reasonable working knowledge of higher 
mathematics would probably take it for granted that a function 
such as F(s;n), which is clearly well-behaved for n %1 and 
Re(s)>0, should possess a Taylor's (or other type of) series repre- 
sentation. Not so. The sleeper in this argument is the condition 
"Re(s)^ 0" . F(s;n) possesses an essential singularity everywhere 
in the left-half complex plane; furthermore, it is not analytic, 
and cannot be continued into the left-half plane. In brief, there 
exist well-behaved functions which can be represented only in 
integral form, and the HGF's encompass a wide class of these, at 
least formally - 
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IV. Derivatives 

Essentially the only simple general formula which can be 
given concerning derivatives of HGF's is 

<i ,n { P V a i' a 2 VV b 2 V x, 3 

(ai> n < a 2V--<*p) n 



• n F_(a,+n,a,+n a n +n * 



(b x ) tb 2 > n ...(b ) - rg-x-'-z™' "P 



b x +n, b 2 +n,,.. f b g +n;x). (4) 

The formulas given below for jFjt^) and 2 F it*' can be gener- 
alized for any particular p,q pair with varying degrees of alge- 
braic effort. Their development involves principally ad hoc 
methods for manipulation of series, and is difficult to systematize. 

(i) For p ■ q = 1, we have 



(a) ( 



£"{*** A **»£ 



= (a) n x a-1 x F y (a+nrb;x), (5) 



(b) tf'rA^i 



{-l) n (b-a) n _ x 
— 2 e x jft (a : b+n;x) , (6) 
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<c> ( 



d Y n f -x b-a+n-1 ) 

■$ ( e x iFi.(«ib,x)J 



- (b-a) n e~ x x ^(a-njbrx) , ( 7) 



(dj (. 



if 1 ! xb_1 1*1<«*>*>J 



(-1) (l-b) n x ^(arb-nTX) . , 8) 



Those four relations demonstrate how either a or b can b< 
increased or decreased by an arbitrary integer, if we wish to 
change both a and b, the requisite operators are concatenated 
(see below, under Continguous Functions ) . 

(ii) for p - 2, q . l, the analogous four relations are 



(a) (*-> 
dx 



n jV*"- 1 2 F lfa ,b;c : xH 



a-1 



ta, n x 3 F (a*n r Hfo>x) , f9 , 



(b, <J|,» [(l-*,» +1 « 2 F l(a . b ,c ;Jt) ] 



(c-a) (c-b) 
~- 2 (l-x) a +»-c-n 2Pi , a , b!C+n , x) _ (10) 



[C > <lf> 



- [««« +n - 1 (l-x,»^« 2 r 1 (a. ) ,,c^ 



(e-a) n x c - a - 1 ,i. It) «*'-c-n j^,^ (u) 
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(d) ^) n ^ c "Vi**' b '°' Jt >y 



■ (c-n) x F (a,b:c-n;x) . (12) 

The last comment of the previous subsection applies here. 



V. Integrals 

Most indefinite integrals which can be evaluated exactly in 
terms of HGF's are simply inversions of the results of the previous 
section. For example, the integral 



r 



x e 9 * p F q (A;Bjx)dx 





cannot in general be represented as a ?-(■)■ 
(i) Two General Integrals 

A vast number o£ definite integrals involving HGF's are 
tabulated {see Bateman (3), Tables of Integral Transforms , Vols. I, 
II). Two typical, and fairly general, examples are 

jV"V rt F (A,B,Xt)dt 
o 

■^ jh-iVV^ V rr b i' b 2 *yA./«>, <13) 

(IAPIACE TRAHSFORM) 
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r l y-1 s-1 

J <l-x) X _F_{A;B?lx)dx 



p q 

= r w+V) p*iVi' a i' a 2 v sib i' b 2' 3+ "^ K (14) 

(MELLIN TRANSFORM) 

(ii) Example: 

We present an example below which demonstrates the power of 
the Hypergeometric formulation. The integral under consideration can 
be evaluatedby more conventional techniques, but only by means of 
"intuitive" reasoning. 

Consider 

I = J x k erf(x)e" ax dx, 
o 



where erf(x) is the Error Function; 

a 2 f* -t 

~WJ_ ■ 



erf (x) * -==• I e ' dt . 



It happens Co be true that (see Tabulation , below) 

2 

Substituting, 



2 f k+1 -(a+l)x £ _ ,„ mMt 2 



I - f±= I x" Ti e v F <l,3/2;x*)dx . 



ftr s Q ll 
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Let x 2 -»t: 

<Tv-J Q 1 1 

But this is identical in form to {13)- Plugging in, we obtain 



Many readers will be singularly unimpressed by this result, inasmuch 
as one problem (the evaluation of the original integral) has been 
replaced by one equally obscure- However, there are several replies 
to this objection, namely; 

(a) The integral has been reduced to a form which is 
easily computable numerically, should the need arise. 

(b) For the particular case, k = 1, one of the members of 
the "A" list is equal to one of the members of the "B" 
list (■ 3/2), so that the F (,) reduces to a iF (.) 
and the result is trivial (see Elementary Properties , 
above) . 

(c) The function 2 F l(*' ^ s subject to a number of variable 
transformations ( Analytic Continuation , below) and 
simple manipulation of the parameters A and h (Con - 
tiguous Functions , below) . These may be invoked for 
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both numerical and theoretical requirements. 

One of these happens to transform the second 
expression on p, 16 into 

i - A . Hi t *»i g ltf . a, 3/ 2,-i) . 

V-rr a . (a+l)K/2 2 a 

It follows that if k is a positive odd integer , 
then a member of the "A"* lis is a negative integer , 
and that the 2*\{-) * s a polynomial of finite degree 
(see Elementary Properties , above) . 

(d) It turns out that even if X is a positive even 
integer there are ati 11 some tricks left in the bag, 
to be described subsequently* 

(e) If k / 0, 1, 2. . . ,n. . . , we can be fairly well assured 
that the representation given in paragraph (c) above 
or the second on p. 16 is the best available. 

For the record, it is not difficult to show that the original 
integral (middle p. 15), using more ad ho c methods,* can be reduced 

to: 

-1 






*by writing erf (x) in its integral representation, switching to polar 
coordinates and performing the r-integral. 
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For arbitrary values of a and k, no compact closed form is 
known for this integral. 

(iii) It is usually desirable that the argument of the HGF 
be the independent variable, x, itself, rather than some f(x). 
If not, then, in general we perform the substitution y » f(x) and 
simplify throughout* If this restriction proves to be impractical 
then it can probably be eased somewhat to allow arguments of the 
form +ax or +ax . 

(iv) The cases in which the integrand involves the product of 
two or more HGF's lies beyond the scope of the present descriptions. 
When only two p F q (-) 's are involved the integral, if it converges in 
any reasonable sense, can be expressed in terms of the aforementioned 
E- and G-Functions. For three of more HGF's in the same integrand 
closed form results are available only in special cases. 

(Zn particular, the last integral given in Bateman (4) , Integral 
Transforms , Vol. II (and if a "God" integral exists, this has to be 
the closest known approximation) is 

a 1( 



/mn / 
M 



p / • G [fix . _ Jdx 

b). h q J K*\r | a ± d s / 



*-l G k + n.l*»/i / |- b l'- b 2 -VV-VV- 



- bl #-b , . . • r-b^fC,, . . .C ,-b '"^fj! 

.-._! (is) 



- 18 - 



which is impressive sort of in the same sense that a 500-pound 
canary is. G (etc.) is Meijer's G-Function, which includes HGF's 
as a subclass*) 

The limited results of this section are intended only to indi- 
cate the potentialities of employing the HGF notation for evaluating 
integrals. However, a surprisingly large number of integrals which 
arise in applications are of the form of either Bq. (13) or Eq. (14) 

In the next section a listing of Integral Representations for 
HGF's is given which widens the present catalogue of integral 
properties. 



VI. Integral Representations 

Simple integral representations of F (.) exist only for 
jFj^.) and jF.t-) for general values of the parameters, excluding 

the elementary cases nF(\(')# i F n'"'' n F t^'^* For ot ^ er values of 
p and q the two integrals of the previous section may be utilized 
iteratively to define any given p Fq(.) in terms of multiple integrals. 

Therefore we emphasize the (single- ) integral representations 
of F (.) and J (•) in this section. 

(i) Integral Representations of F (a*b*x): 
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(a) jp^aibrx) - 



r(a>P(b 



-; f^iui^.V lie: 



(b) 



, r'(-:-o x i-b/* e v v -i 



to) 



•(x-v) b_a - 1 dv. 



(b-a)P(a) 



^ ."« w b " a - 2 



(17) 



(1-w) 3 " 1 dw. 



(18) 



(d) 



r(a)P(b-a) -ii 



a-s- 



a- 1 



- (i=S) ds, 

1+a 



(19) 



And, in general. 



(•) 



Tlbjexp {- -S*l 



ra>-a)r<a) 



*ds (< u - c) * x (*- u ) b " a ~ 1 



where, in all cases, we assume that Re{b) > Re{a) > 0. 



(ii) Integral Representation of jfi (a,b;c;x) : 
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(a) jj F lt a ' b?c ' x ' 



Cm r >- 



P(b)r*(c-b) 



(l+s) a - c (l+5x)- a da , (21) 



(b) ^(a.b.-crl-x) 



(c) 2 P 1 (a,b;c:l/x) 



no 


*s: 


t b-l . 


■ (i-t) c 


-b" 1 (1. 


-xt)" a dt , 


■*J& 


.. f 


Vi.*-"- 1 



(22; 



■ s* C (s-z) _a ds ( (23) 



(d) 2 F,<a,bjcjx) = rtc) r 

2 i r(br f (c-bi J„ 



P(b)P(c-b) J Q 



_t c-b-1 ,. _- 

" (1-e ) (l-ze _t ) a dt . (24) 



Since a and b are equivalent parameters, another set of four relation- 
can be generated simply by interchanging a and b everywhere. As 
above, we require that Re(c)> Re(b) (or Re(a))>0. 

(iii) There are a large number of alternative Integral Repre- 
sentations, of course, which are obtained by transformation of the 
variables of integration. These have been omitted for the moBt part 
for reasons given in the subsequent remarX. Also, any -,*"„(•) can *>e 
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represented in terms of a complex contour integral known as a 
Mellin-Bames line integral (Slater (1), (2)). These representa- 
tions have also been disregarded, at least for tne present, because 
of the difficulty of establishing simple rules for the selection of 
the proper contour of integration, as well as the lack of a mechan- 
ism for performing Residue Calculus. 

(iv) Throughout this Memo it has been implicitly assumed that 
the independent variables of the BGF are simple functions. When 
this is not the case, especially in applications involving integrals 
and differential equations, it would seem to be more efficient to 
first make an appropriate change of variables and then apply a rela- 
tively small amount of sorting to the result, as opposed to having 
to deal with a much longer list of more general forms. 

For example, given 



C^n , . _2b-l , _2c-2b-l 
! = j [si* t) (cost) dt ( 

^o (1 - z sin 2 t) 



we set x = sin (t) and obtain after simplification: 



r 



. I x*" 1 (l-x) C - b_1 <X-zx)~ a dx, 



■ 2 Hfflj-b) 2Fl (a, biC;Z ) , 
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according to formula (ii) (b) of this section. 

This reduction is not trivial to recognize, of course. The 
motivation is to seek out the "most troublesome*' term in the 
expression — the denominator of the integrand — and break it down 
into simpler terms, in the hope that the procedure does not make 
an even worse mess of the remaining elements of the expression. In 
this sense, it is difficult to define "most troublesome" in a succinct 
manner. For this particular case, an equivalent but less direct path 
for the reduction would be to let 

y o sin(t), then x ■ y , or 

y = coe(t), then x = y 2 , 
assuming we "know" that sin 2 (t) ■ 1 - cos z (t). 

On the other hand, given an integrand which involves only trigo- 
nometric (or hyperbolic) functions of the independent variable, the 
substitution y =* sin(t), or cos(t) (sinh(t), or cosh(t)) reduces it 
to a purely algebraic expression; if further reductions are required, 
they shotild then be relatively easy to spot. 

VII. Differential Equations 

Tho discussion in this section has been purposely abridged 
principally because an exhaustive cataloguing of the homogeneous 
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solutions of DE's of the Hyper geometric type would tax the patience 
of the reader. Furthermore, these are well documented in References 
(1), (2), and (3), Therefore, the remarks below are principally of 
a general nature. However, efforts bb being made at the present 
time to reduce this detailed classification to a less unwieldy set 
of fundamental procedures and operations. 

The most general (and essentially only and basic) solution 
technique for linear DE's with polynomial coefficients is that of 
Frobenius (Slater (1), (2)). Other, more concise methods are 
nearly always limited in scope and are generally equivalent to 
Frobenius ' in any case. 

The advantage of the HGF formulation is the fact that the behavior 
of the solutions with respect to the nature of the singular points 
and other types of limit properties are compactly "encoded" into the 
"A" and H B" lists. 

(i) The General Hypergeometric Differential Equation: 
The series D F a (A:IJ;x) satisfies 



[■ 



-^(x-£ + b. - 1) <x-* + b 2 - 1) ... (x-1 + b_ - 1) 
dx dx * dx * dx q 



- x(x-^ + a L ) (x-£ + a ? ) ... (x-1 + *J\ y - . (25) 



dx x dx 2 dx P 
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The order of this DE is max{p,q+l). 

For further elaboration and classification of singular points 
see Slater (2) , pp. 42-45. 

(ii) Gauss ' Equation: 
The aeries -> F \ ( a * k? c; x) , known as Gauss' HGF, satisfies 

dx : 



x(I-x)4Jf + T c _ (i +a +b)x"l |y - aby > . 
dx 2 L J ax 



(26) 



See Slater (2) , pp. 5-13. 

(iii) Kurcmer's Equation: 
The aeries jPit.) (Rummer's HGF) satisfies 

x^ + (b-x)ft _ ay = o . {27) 

dx 2 ■■ 

See Slater (1), Ch. 1. 

(iv) A Class of Second-Order Linear DE's: 
Suppose we are required to find the solutions of the homogeneous 
equation 

P 2 (x) ^¥ + P,(x) |* + Po (x)y => , (28) 

dx 

where P2 (x) , Pi(x), and Pq(x) are real polynomials in x, 
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We inquire: can the solution(s) we written in terms of, say, 
g(x) '2Fita,b;cjpx + q), 

where g(x) is an elementary function? The solutions might also 

involve 2 F 0**'' 1 P 1* - '' or F 1*"'" Then a 9 ain * ifc ^V be true that 
the DE is not of Hypergeonetric form. The details for transforming 

any given DE into the latter have not yet been completely worked 
out — in fact, to the best of this writer's knowledge, it has never 
been attempted — and there is a strong possibility that it may be 
sijnpler in general to automate Frobenius ' method- On the other 
hand, the great majority of DE's of this form which are encountered 
in applications possess solutions in terms of 2 F l(~) or i^i(-)- 
By testing for the location and nature of the singular points of the 
given DE, which is a fairly straightforward task consisting mainly 
of finding the roots of the p^(x), the possibilities can be cut down 
significantly also. 

{v) General Linoar DE with Polynomial Coefficients: 
A great deal of additional work is required in the area of the 
investigation of HGF's in trie DE context. This is not to say that 
most of the results are not known — rather, they are numerous and 
widely dispersed throughout the literature and require systematization. 
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The best efforts to date are due to Slater {(1) and (2), but there 
is much left to be desired, especially where the aims of MATHLAB 
are concerned. 

Suppose wo are given 



tiXii) 



Y - , (29) 



where the p, (x) are real polynomials in x. The conditions under 
which this DE can be transformed into the Hypergeometric Equation 
are not known in general for n * 4, and only in principle for the 
cases of n - 4 or 3. (The latter require the solution of systsns 
of algebraic equations of degree 4 and 3 respectively.) As before, 
by locating and classifying the singularities of the equation, it 
may be possible to determine if the given DE is not of Hypergeometri 
form; however, if we cannot prove that it is not Hypergeometric, it 
may still be or not be Hypergeometric- The Hypergeometric DE pos- 
sesses a£ most three distinct singularities: if the independent 
variable is x these lie at 0, 1, and infinity. The singularity at 
the origin may be reular or irregular; that at 1 is regular or 
removable; that at infinity is regular or irregular* 
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VIII* Analytic Continuation 

Analytic Continuation Formulas are known completely only for 

the HGF's JV. (9)4 l'lt*)* an ^ 'he trivial case iFq(.)* For other 

F (O's these can nearly always be generated in any particular 
p q 

instance, but only by resort to fairly sophisticated techniques 
which require a thorough background of complex variable theory, 
and thus lie beyond the limited scope of the presont study* 

(Analytic Continuation and the closely related topic of Con- 
tiguous Functions, discussed in the following section, are useful 

principally for numerical applications and also for the investiga- 
tion of asymptotic properties. See remark (iv) below.) 

(i) Trivial Case: 

,F ft <a; ;x) - (l-x)~ a - e" i1la x" a (1 - lA)-* 

.*-(*; ;l/x) . (30) 



l r 0' 



x* 



r 



(ii) Kummor's Theorem: 

jFjfajbix) = e x .Fitb-ajb.'-x), or 

1 

e x ^(arb.'-x) = ^(b-aibrx) . 
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(iii) Gauas-Euler Reductions for Gauss' HGFi 

(a) ^(a.b^x) = (l-x) c --b ^(c-a.-c-b^rx) . 



!32) 



(b) = (1_ X ) 



2 F 1 ^,C-h; X: -S.) , (33) 



+ (l-x) c - a - b C (OH a+b-cl . 

Rairfb) 

• 2 F 1 (c-o,c-b f c-a-b+l;l-x) , (34j 



(d> ^ "rtbMcIaJ <-*>" a ^(a.l-c+a, 1-b+a • 1/x) 

* ruif%2} ( - x) " b A< b -***» i-«* vx, , 



35) 



,oi •■ (1 - c '- fM^f ^ t - c - bt - b * i 'A' 



X 

(36) 
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» ■ R£fei *- ^l-.*-W*HfcO*,l - i, 



• ^(c-aa-afc-a-b+lrl - -i) - (37) 



The conditions under which those transformations are valid is 
usually equivalent to the requirement that the arguments on both sides 
of the equality are defined, and that the various Gamma Functions do 
not possess zero or negative integer arguments. 

(iv) General Remarks: 

Analytic continuation is applicable only when p = q+l f or p = q. 
For the case where p> q+1, the HGF is not analytic, and when p<q the 
series converges uniformly everywhere (see Elementary properties (i) 
above) - 

The usefulness of Analytic Continuation arises in both numerical 
and theoretical areas* For example, we may wish to evaluate a 

F-(a,brc;x) for large values of x- Direct substitution into the 
series may lead to an unwieldy calculation, due to slow convergence. 

However, we can transform: x >-l/x {see previous article) and the 

resulting series will converge relatively quickly as x— h« - On the 
other hand, the purely analytical advantages are well demonstrated in 
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the example given above (V(ii)) where the tFj (.) which emerged from 
the woodwork {line 4 p. 16) was transformed into a more readily inter- 
pretable form (line 4, p. 17). 

Another valuable numerical and theoretical tool arises frccn the 

contiguous properties of RQF'Sj discusscc in the next; soction. 

IX. Contiguous Functions 

The two HGF's p F (A'?B'tx), F (A" ;B" jx) are said to be contiguous 
when all parameters contained in their respective "A" and "B" lists 
are equal, except for one pair of parameters (from either A or B) , and 
these two differ by unity . 

For example, 

2 F ,a l' a 2' '*>' 2 F 0< a l' a 2 +1? JX * 
1 P 2 (a 1 ;b 1 ,b 2 fx) ^(ft^b^-l.bgix) 

are contiguous pairs, but the pair 

is not. 

However, the last pair of functions are said to be associated . 
In general, two -JPL (■) 's a** associated if the respective elements 
of their parameter lists (both A and B, or either) differ by integers 
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or zero . 

It is not difficult to show from the differential properties 
(section IV) that a linear relationship exists between ro+1 contiguous 
P 's, where m is equal to the order of the associated DE {Eq. (25)), 
i.e. ■ max (p,g+l). For example, for 2 Fi(-) or iFi(0* the number of 
functions in the expression is three . Furthermore, it can be shown 
that a similar three- function relation exists between any three 

associated 2 F l' s or 1 F 1 S * 

The term "linear" in this context means that the coefficients of 

the various p F q 's in the relationship are linear functions of any 

specifically chosen parameter, or the independent variable. For 

example, one (of six) contiguous identities for F.{-) is 

b • jF^arbyx) = b - jFita-lrbjx) - x • jFiUrb+l-x) . 

The implication of this result is that, in the cases of 2 F lt*J 
and iF^(-) at least, we can compute: 

(a) ? F, {a+m, b+nj c+1; x) from 

2 F i (a,brc;x) and 2 F l* a ' b:C+1?x ) * 
for example. 

(b) 1P1 (a+mrb+n;x) from 
2 F 1 (ajb;x) and jF^ajb+lrx) , 
for example^ 
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where m, n and 1 are integers, positive or negative. 

For other values of p and q, it is generally simpler to iter- 
ate co ntiguous relationships in order to calculate an associated 
function. 

Complete listings of the contiguous relationships have been pre- 
pared only for -F.(.) and iF*(*-)« These are not given here for rea- 
sons of space, but may be found respectively in Slater (2), pp. 13-14, 
and Slater (1), p. 19. 

It should be remarked that the derivation of contiguous relation- 
ships (i.e. the determination of the coefficients) for general values 
of p and g is by no means trivial. The calculation requires fairly 
careful handling of power series relationships, and there is apparently 
no general systematic technique. 

We recall for a moment the remark made above. Section V (ii)(d), 

concerning "tricks": all o^i >& °^ tllG form appearing on that page, 
where k is a pos it i ve even inteqer, are associated functions. It 
follows from the results just presented that if we know the value of 
the integral for k = and 2, say, we can generate in a simple manner 
the values for k ■ 4, 6, 8, ... 2n, ■*•* In fact, this starting pair 
is easily calculated frcro the integral on the bottom of p. 17 by ele- 
mentary methods. 
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X. Tabulation of Functions Representable in Hypergeometric Form 

In this section we list some of the more widely known Higher 

Transcendental Functions which can be expressed aa F {.)• 

P *J 

The trivial cases ,F-(.) and F (.) have already been mentioned 
(Section III) and will not be repeated. For completeness, we should 
include in this category the rather obvious identity 

pFptA;Arx) = e* . 

The reader is cautioned to be wary of notation; the one employed 
here is generally that of Ref. (4), since it seems to represent a 
reasonable compromise among the {articular variations which are encoun- 
tered in a cross-section of the different scientific disciplines. Also, 
Ref. (4) nearly always catalogues the most common variant notations for 
each of the functions as they are introduced - 

(i) Functions Ropresentable by iFi(a;b;x)r 

(a) Bessel Functions 

r (ay2a 7 x) - Pl^e^ 2 (x/4)^ a i^(x/2) . (38) 

Depending upon whether a is an integer or half of an odd integer, 
positive or negative, and whether x is positive, negative or imagi- 
nary, a variety of differently named members of the family of Bessel 
Functions can be generated. (See Ref- (4), Eqs. 13.6.1 - 13.6.7.) 
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(b) Polynomials 

jFj^-njbjx) - -£l- L n (x) - (LAGUERRE) (39) 

When b = 1/2 and 3/2, we obtain respectively the even- and odd-order 
Hermits Polynomials: 



l'l«-»'^.X) - ^T t-W- H- 2n+1 i/H, .j ,40, 

(BERMITE) 
Another type, essentially the same as (39), is 

^(-nr-^l-x) B -|^I__^_p n (V, X ) . (41) 

(POISSON-CHARLIER) 



(c) Miscellaneous Functions 

jF^ara+lj-x) - ax~ a if (a,x) . (42 ) 

(INCOMPLETE GAMMA) 
When a = 4, this reduces to the error function! 

jE^ttfS^r-x 2 ) -*nTerf<x) . (ERROR) (43) 

A generalization of the Hermite Polynomials, just cited (40), is 
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!?!< -*T/2,1/2 JX ) =^e x / 2 E v <°> (/2x) . 

1 F 1 (V3 - *^2; 3/2; x) = -^ e 1 ^ 2 E^ 1 ' (J25) . 

(PARABOLIC CYLINDER) 
F (H+i. n+ l : r2) - n:r- 2n - m -^ e r 2 T{m , n ,r) . (45) 



fW 



(TORONTO) 



It diould be borne in mind that Hummer's Transformation (£q.(31)) 
can be applied throughout, essentially doubling the size of the list. 

(ii) Functions Repreaentable by 2 F i (a,b;c;x) : 
(a) Polynomials 

(JACOB I) 
(Clearly, any 2 F ] (-) which reduces to a finite polynomial can be alter- 
natively represented by a Jacobi Polynomial by means of a proper choice 
of n, *r, and J. Furthermore, by employing the Analytic Continuation 
formulas of VIII (iii) at least five more representations equivalent to 
(46) can be obtained,) 

C n ( " W - £^ 2 F 1 (-n.n.2«, **, if) . <47, 

(GEGENBAUER) 
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1-x 

t 



r p (x) - 2 T x {-n,n-, hi ±=2- ) 

(48) 



D n (x) ■= (n + l). 2 Fi(-n,n+l ; 3/2, ±z* ) 



(TSCHEBYCHEV) 



P n (x) = ^l-n.n+l; 1, I|£) . (LEGENDRE) (49) 



The results of (47)- (49) are, of course, subject to the transfor- 
mations of VIII (iii). where valid. 

(b) Elementary Functions: 



2 P 1 {l.ll2;x) - log (1-x) 



2 F 



^ftrifS/a*-**) -pt* (X) • 



t- L"i l-2a 



(50) 



lfcl|V *rf> . ^ log(^). ("» 



-P. A,»»r3/2»X 2 ) --isin" X (x) . (52) 



-^an^tx) . <53) 

«i,-2a „ >h-3a 



(54) 



(55) 
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F. (2a,a+l;a;x) = 



1+z 



2'1 



(1-z) 



7STT 



(56) 



(Notice that the last three results can also be expressed in terms 
of simple functions of jFQf.J's-) 

(c) Other Functions: 






(57) 



(COMPLETE ELLIPTIC INTERVALS) 



"P 



^(P.l-qrp+ljx) - px r B x (p,q) . 

(INCOMPLETE BETA) 



(58) 



,F, (a,b;c;x) - — I tg) — E(a,b;c;--t . 



{MACBOBERT'S E-) 



(59) 



r<a)r(b> 22 ^ x a>b y 



(60) 



(MEITER'S G-) 
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(1 + x) n to N terms: 



- 1 + (J)x+ fyx 2 + ... + (JJ )X « 



i F o'^ n? ; ~ x ' to N ternifi 



{TRUNCATED BINOMIAL SERIES) 
(See generalization in next section . > 



(iii) other Representations: 

For p F q* s other th * n the cases already listed, the reader ie 
referred to Slater (2), pp. 46-7, and Bateman (4), Higher Transcen- 
dental Functions, Vol, II. under "Generalized Hypergecwnetric Series" 

One result which is of interest, much like that of Bq- (61) 
above, is the following: 

« x to N terms 



" Q F o ( ; ;xJ to N t 



ernis 



x N 1 

N . 2 P <-»»!■• •- ~> ■ (62) 

(TRUNCATED EXPONENTIAL SERIES) 
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(Notice that although it was stated previously that if p > q + l 
the associated D F a (.) is undefined in general, this restriction 
does not apply if one or more members of the "A" list are nonposi* 
tive integers, in which case the series terminates in a finite 
number of terms, as in the present example-) 

All other truncated series of circular and hyperbolic func- 
tions can be represented in a like manner. For example, 

cosh x to N terms 

= H(* x + e" X ) to N terms 






XI. Miscellaneous Results and Random Eacts 

This section lists some useful and powerful properties which 
do not fit neatly into any of the above categorizations. 

(i) Truncation of Series 

Suppose that we are given a power series whose individual 
terms are of the form of those of a p F q (*)# but which is terminated 
after K terms, as in Eqs. (61) and (62) above. We state the follow- 
ing theorem: 
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[pV* ; £ ;x) ]n 



N 



f Y2 < a l>k< a 2>k---( a p>i c x* 
k=0 ( fa l>k (b 2>k'"* b p*k * : 

l_ H _A f H>*»* . , (63, 



where 1-n-B*-» [1-B-bj, , (l-H-b2) <l-N-b q ), 

1-H-A^-Kl-N-ai, , (l-N-a 2 ) (l-N-a p ) . 

Eq. (63) is a generalization of the aforementioned examples. 

(ii) Reversal of Series 

A closely related result applies to terminating Q F 's, i.e. 
those which possess or a negative integer in their "A" list. The 
aeries can be expressed in terms of powers of (1/x) as follows: 



p P q {a 1 .a 2 ,...a p _ 1 ,-m f bi,b a , . ■ .b^j x) 

-(b l)m (b 2 , m ... (bq P , m <"*> * ^ F M'H. !-»» 



i_a; t=agsa , . {64 , 
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(iii) Particular Values of the Argument 

For certain F n ' 8 > especially when p m q+l f the series is 
summable in closed form when the independent variable takes on 
special values- some of these are listed below. 






3 P,(a,br 1+a-b; -1) - 2" a jJl+a-b)T(l * V 2 ) 



(65) 



/2) 



(66) 



2 F 1 (a,l-a; c; h) - — ■ g . (67) 

2 1 p ( £±£jp(i±£zl) 



A few others for F (.) can be generated by means of the 

Analytic Continuation formulas given above. For other pFq' s °^ ***e 

form „F (»)- with particular arguments, see Slater (2), pp. 48-84. 
p+1 p 

(The underlying theory of the latter is quite deep, and is of purely 
abstract interest for the most part, in any case.) 



(iv) Confluence 

The adjective "confluent" is often encountered in the literature 
of HyjE rgeametric Functions. The word means, literally, "running 
together". In the present context, the term denotes the "running 
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together" of the singular points of a F (•) at one and infinity* 

p q 

To fully appreciate the manner in which this is carried off would 
require a lengthy and complicated discussion of the DE theory of 
HGF's. For example, if we make the change of variables x — *x/b 
in Gauss' Equation (26), and then let b — ±m (carefully) we obtain 
Rummer's Equation (27) (with b replaced by c) . 

Roughly speaking, 2 F l* a ' b ' c ? ar 9) possesses a regular singu- 
larity when arg « 1, as well as the singularities at arg ■ and 
<»; thus, if arg ■ x/b, the "middle" singular point lies at b, 
and thus goes to infinity as b-*aj (to join the one already there). 

All of this is easier to visualize in terms of power series; 
suppose we want to "do a confluence" on a given o Fq(0 with respect 
to a member of the "A" list, say a ( the last membor) . We can 
represent the operation symbolically as 

CONFLUENCE (pF q (A;Brx)) wrt a 

x 



- Lim (pF q (a 1 ,a 2 , •■•*]/*' £; — ) ) 



■uv*- p - X P a P 



' p-l F q (a l' a 2'-"*p-l ; & x) * < 68 > 



This result is fairly obvious if one recalls Eq* (1): the 
only terms in the summand which involve a p constitute the ratio: 
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"V" V fl p »• a P 



1{1 + l/a p )(l + 2/a )...(! + »=i) , 

a P 



which clearly goes to unity as a_— +*w, thus wiping out all trace 
M a and leaving the independent variable unscathed* 

Thus, briefly, to obtain a confluent form of a F (,)■ one 
member of the "A" list is deleted and p is reduced by one. 

By way of justification of all of this verbiage, which may on 
the surface appear to be describing a fairly trivial property, let 
it be noted that Slater (1), (2) employs confluence nearly every- 
where to demonstrate theorems regarding jF^'s by "doing a confluence 1 
throughout the analogous theorems for 2 Fi- (The overall theory of 
the latter is richer and more versatile.) in the same vein, but in 
reverse, this writer has on several occasions dealt with a recalci- 
trant jFjtO (say, in an integrand) as follows: 
Let 

F (a;c;x)— > jF^a^-crx/h) , 

play 2 F l" type sanies, get an answer (if possible), and then after all 
the dust has settled let b— *<& . if the result makes sense at all, 
it is probably correct* 
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(v) Remarks 

There area large number of other isolated properties of HGF's 
which are not described here in detail for reasons of space and 
present applicability. The most important of these concern: 
Addition Theorems , i.e. the expansion of D F_(A;Ji;x+y} in terms of 
pFq's involving x or y singly; similarly. Multiplication Theorem a 
for F (A;I3;xy)r expansions of D F q (,) in terms of HGF's with smaller 
values of p and/or g. (As an elementary example of the latter, 
recall that X*r(x) can be expressed alternatively as -FJ.) or 

For further documentation, the reader is referred to slater (1), 
(2), Bateman (3), or the Bureau of Standards Table (4). 

XII. _ Important Omissions 

At several junctures in this report this writer has passed 
over or detoured subject matter which should properly have been 
included, but for the fact that the text would have at least doubled 
in size. Considerable pains have been taken to sift through the body 
of theory in order to simultaneously maintain readability, pertinence 
and brevity, and wherever gaps or incomplete listings have arisen, a 
specific reference has generally been provided. 

However, at least two significant topics have received inadequate 
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attention. 

(i) Differential Equtions: 

As indicated in the closing remarks of Section VII above, 
the theory of RGF's as applied to linear ordinary DE's with poly- 
nomial coefficients requires extensive study- A tabulation of 
the known results was not included because, at least to the insight 
of this author, no "pattern" seemed to emerge, except at the elemen- 
tary "Frobenius Method" level, 

(ii) Associated Kummer Functions: 

Associated Kummer Functions constitute a class of solutions 
of Kummer's Equation (27) for which "i^il-)" is not valid- in 
many cases, these functions {usually denoted by U{a;brx)) can be 
written in terms of .P, 's? when they cannot, they are nearly always 
characterized by logarithmic singularities or other types of branch 
points near the origin. 

The principal reason that these were not mentioned is that 
the associated theory parallels that of ,F*(.) except for a number 
of nagging exceptions which are just sufficient to create confusion 
{very much like attempting to match up trig- function and hyperbolic 
function identities*)- Eventually, of course, these functions will 



* cosh 2 x + sinh x - ? (quickly, now) 
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have to be included, but it was felt in the course of the prepara- 
tion of this report that the penalty in terms of length and reader 
attention did not justify the relatively minor virtue of a more 
complete presentation. 

A listing of functions representable by U(a;b;x) is given in 
Bef. (4), p- 510. It can be seen by comparison to the list on 
p. 509 (for ^j^arbjx)) that the degree of nonoverlap is not serious. 

XIII. LISP Representation 

R. Fateman has suggested that the following format be adopted 
for LISP representation of hypergeometric functions: 

F (A?B;X) 
p q 

«— | ((MFHYP P 0)(A 1 A 2 ... OtBi * 2 * " B q ,Xj (69) 

The car of the right side identifies the function as being a 
F (.); the cadr and caddr are the "A" and "B" lists, respectively, 

p q 

and the cadddr is the independent variable. To be slightly more 
general, the last element should probably be replaced by ("FUNCTION x) 
or a quoted lambda expression, in the event that the independent 
variable is not x itself. 

The inclusion of the parameters P and in the car is admit- 
tedly redundant, inasmuch as their values are specified by the 
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length* of the "A" and "B" lists; however, it appears at this time 
that it is more convenient in the long run to have P and Q 
available at a fairly high level in the form of atoms rather than 
cluttering up the landscape with expressions like (LENGTH (CADR Y) > 
and so forth. Furthermore, the decision paths which determine how 
a particular F {.) is to be manipulated branch at a very early 
stage, depending upon the values of P and 0, somewhat as follows. 

After reduction (i.e. cancellation of common members of the 
"A" and "B" lists), we examine P and Q to see if Pq(*)- q F i*'* 
or F (.) has emerged; if not, then test for jFiU): if T, then 
branch to a recognition package for jF^C*)' if not < t#8t for 2 F l'-' ; 
if T, then branch to a recognition package for 2 F i t *' r otherwtsei 

(a) If "A" lilt contains or negative integers, and "B" list 
does not, the F (.) is a finite polynomial. 

P H 

(b) If "B" list contains or negative integers, and "A" list 
does not, function is undefined. 

(c) If both "A" and "B" lists contain or negative integers, 
and all such members of the "A" list are greater than all such mem- 
bers of the "B" list, then function is a polynomial; otherwise 

undefined. 

(d) Otherwise (i.e. none of "A M and/or "B" lists are or 
negative integers) function is defined in terms of an infinite 
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series if p^q+1, undefined otherwise. 

A flow chart for this logic is shown in Pig. 1, by way of 
example. Other packages for integrals, derivatives, and differen- 
tial equations have also been flowcharted. (See remarks in 
Conclusions , below-) 

The overall input is assumed to be a HGF which has emerged 
from lowest-level operations in "unrefined* 1 form, i.e. the package 
is presented with a F (.) with no a pri ori knowledge. We call 
this input G, and the output of the reduction routine F. The nota- 
tions "T" and "NIL" on the branches denote M yes" and "no"; a, and 
b, denote members of the "A" and "B" lists. 

(i) HYPREDUCE 




Fig, la . (G* denotes the reduced form of G which omits a* 
and b., if they are equal and reduces both p and q by 
unity.) 
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(ii) HYPSORT 



r N^ 



NIL 

p ■ 2\ j HYPSORT3 



HYPSORT0 



HYPSORT1 



'I 



HVPSORT2 



Fig, lb. 



!iii) HYPSORT0 




NIL 

2J — *E 



_ NIL. 

3— 4 



LCGKEZSSEL 




Fig- lc . LOOKBESSEL looks at F (which is, by now, a P.(.) f 
and therefore a Bessel Function) and determines which°parti- 
cular Bessel Function best fits the parameters. 



NONPOSBP (prodicate) examines the "B" list; T is returned if 
none of the membersareO or a negative integer? otherwise NIL, 
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(iv) HVPS0RT1 



NIL 



NIL 




Fi 9- ld - LOOK1F1 i8 a recognition package for 1F1 (.) which 
attempts to identify F in a list of known representations 



(see X(i), above) . 



(V) HYP50RT2 



Hi, 



ni: 



HYP SORT 3 



LO0K2F1 



Fig- le . L00K2F1: like LOOKlFl, but for _F. ( . ) (see 
X(ii) , above) . 2 1 
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(vi) HYPS0RT3 



I ^ | p > q+l| 




Fig, If , NONPOSAP (predicate) is identical to NONPOSBP, 
except operates on "A" list . 

FLANKABP (predicate) flags members of both B A" and "B" lists 
which are or a negative integer: returns T if all such mem- 
bers of M A" list are greater than all such members of "B" 
list: otherwise NIL. 



Some fragments of this package and of others have been coded 
and appear to work satisfactorily* The most difficult problems in 
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the programming area are anticipated to lie in the interfacing 
with the present MATHLAB, particularly in the "simplification" 
procedures, and in designing conventions for I/O which are effi- 
cient frcra the point of view of both the user and the system. 

Host of these problems have not yet been considered in detail, 
although none of them appear insurmountable at present* At the 
moment, other more urgent problems are at hand, so that this work 
will have to be suspended temporarily. 

XIV. Conclusions 

The purpose of this report has been to introduce the reader to 
Hypergeometric Functions as painlessly as possible, and to give 
scroe indication of the strengths and weaknesses of the Hypergeotnetric 
Formulation. 

Its power arises from its generality and the fact that it 
enables one to deal with a much wider class of functions than can 
be handled at present, and to do so in a unified manner. The several 
more-or-less detailed examples which are scattered throughout the 
text were chosen to support this contention. 

On the other hand, implementing these techniques presents a 
large number of systems and programming problems, and it would be 
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difficult at tho present time to predict how much time and effort 
would be required to solve them. Furthermore, such a package would 
probably be costly in terms of core and running tijne, but that is 
of secondary importance- 
It is this writer's firm belief that if MATHIAB is in the 
future to appeal to a wide variety of users, fran freshman calculus 
student to theoretical physicist, it will have to eventually incor- 
porate some or all of the methods outlined in this report. There 
is literally no other sensible way of dealing with Higher Trans- 
cendental Functions efficiently within the confines of a finite 
system. 

Perhaps this study is premature, inasmuch as there are bugs 
in the present system which must be eliminated, as well as some 
new foundation blocks to be incorporated before attempting to build 
in the Hypergeometric Function methods. At least, these have now 
been documented for future reference and a promising direction for 
research has been provided. 
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